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We  have  developed  computer  program  which  calculates  the  sea-surface  three- 
dimensional  gravity  effect  of  a conical  seamount  from  a point  above  the  apex 
to  any  distance  off-axis.  The  results  of  this  program  are  compared  with 
those  of  a Taiwan!  two-dimensional  equivalent  of  a ship  track  across  the  off- 
axis  hyperbolic  slice  of  the  cone.  A study  of  more  than  100  seamounts  in  the 
Pacific  Ocean  showed  that  their  slope  angles  range  from4.5*  to  17  . For  this 
reason,  model  cases  with  slope  angles  of  5*,  10*  and  15  in  water  depths 
of  1 to  6 km  were  considered.  For  a seamount  with  its  apex  at  2-km  depth 
and  base  at  5-kn  depth  and  10*  slope  angle,  the  two-dimensional  method  over- 
corrects by  18 % over  the  apex,  but  gives  the  same  result  as  the  three- 
dimensional  method  at  an  off-axis  distance  of  60£  of  the  base  radius.  Farther 
out  than  60t  of  the  base  radius,  the  two-dimensional  method  under corrects. 

We  have  developed  a second  computer  program  which  uses  a two-dimensional 
method  to  correct  the  free-air  anomaly  for  topography  along  the  ship's  track. 
Curves  from  the  cone  program  as  well  as  two-dimensional  topography  corrected 
free- air  anomaly  tracks  are  presented. 
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ABSTRACT 


Localized  bathymetric  features  produce  large  amplitude 
perturbations  in  the  marine  free-air  gravity  anomalies.  In  order 
to  study  long  wavelength  anomalies,  which  are  of  particular 
interest  in -marine  geodesy,  it  is  desirable  to  remove  the  topo- 
graphy "noise  level".  However,  to  do  this  in  a truly  three- 
dimensional  fashion  along  a marine  gravity  track  is  not  practical 
because  of  the  formidable  amount  of  data  preparation  required. 

It  is  also  undesirable  to  use  conventional  Bouguer  anomalies 
because  these  can  easily  he  in  error  by  607.  because  of  the  one- 
dimensional nature  of  the  correction.  We  have  developed  a 
computer  program  which  calculates  the  sea-surface  three- 
dimensional  gravity  effect  of  a conical  seamount  from  a point 
above  the  apex  to  any  distance  off-axis.  The  results  of  this  . 

program  are  compared  with  those  of  a Talwani  two-dimensional 
equivalent  of  a ship  track  across  the  off-axis  hyperbolic  slice 
of  the  cone.  A study  of  more  than  ICO  seamounts  in  the  Pacific 
Ocean  showed  that  their  slope  angles  range  from  4.5*  to  17*.  For  . 
this  reason,  model  cases  with  slope  angles  of  5,  10  and  15  - 

in  water  depths  of  1 to  6 km  were  considered.  For  a seamount 
with  its  apex  at  2 km  depth  and  base  at  5 km  depth  a 10*  slop* 
angle,  the  two-dimensional  method  over-correc ts  by  18%  over  th* 
apex,  but  gives  the  same  result  as  the  three-dimensional  method 
at  an  off-axis  distance  of  60%  of  the  base  radius.  Farther  out 
than  607.  of  the  base  radius,  the  two-dimensional  method  under- 
corrects. We  have  developed  a second  computer  program  which 
uses  a two-dimensional  method  to  correct  the  free-air  anomaly 
for  topography  along  thfe  ship's  track.  Curves  from  the  con* 
program  as  well  as  two-dimensional  topography 'corrected  free-air 
anomaly  tracks  are  presented. 


INTRODUCTION 

Gravity  anomaly  data,  particularly  in  the  form  of  l*xl* 
means,  are  needed  in  order  to  determine  undulutions  of  the 
gravimetric  geoid,  determine  components  of  the  deflection  of  the 
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a large  number  of  stations.  Any  station  location  Is  well  know, 
and  the  gravimeter  is  stationary  during  the  observation.  Many 
areas  of  large  extent  have  observations  on  a grid-like  pattern, 
and  there  are  enough  point  anomalies  available  so  that  simple 
averaging  can  be  done  and  a reliable  mean  is  obtained.  Good 
topographic  maps  are  available  so  that  terrain  corrections  can  be 
made. 

The  situation  is  not  so  good  as  regards  marine  gravity 
results.  One  factor  is  the  high  cost  of  the  sea  gravimeter 
Itself,  at  least  15  times  the  cost  of  land  meter.  Another  is  the 
very  high  cost  of  ship  operation.  A ship  suitable  for  sea- 
gravity  observations  costs  thousands  of  dollars  a day  to  operate. 
These  two  factors  alone  indicate  not  only  why  there  are  great 
gaps  in  the  sea  gravity  data,  but  also  why  the  areal  coverage 
usually  consists  of  a long  single  track  through  a region  rather 
than  a grid  pattern  of  observations.  Then  it  becomes  much  more 
difficult  with  the  sea  data  to  determine  mean  l*xl*  anomalies 
than  it  has  been  for  the  land  data.  To  make  matters  worse,  the 
ocean  bottom  very  often  has  relatively  high  relief  which  usually 
consists  of  seamounts.  Because  of  navigation  uncertainties,  the 
exact  location  as  well  as  configuration  of  bottom  topography  is 
not  well  know.  In  addition,  slope  corrections  to  the  echo 
sounder  depths  have  probably  not  been  made,  so  that  in  bottom 
topography  configurations  where  there  is  a high  slope  angle,  the 
charted  depths  may  be  incorrect.  In  this  regard  see  Krause 
(1962) . 

For  land  gravity  observations,  a several  km  high  topographic 
fea  ture  only  a few  kilometers  away  has  virtually  no  terrain 
effect  on  the  vertical  component  of  gravity.  At  sea,  however, 
the  observation  point,  instead  of  being  at  almost  a right  angle 
to  the  feature,  is  usually  above  the  level  of  the  top  of  the 
fea  tore.  Therefore  the  angle  from  the  observation  point  to  the 
center  of  mass  of  the  topographic  feature  will  be  considerably 
different  from  nearly  horizontal,  and  the  gravity  effect  much 
larger. 

It  is  well  known  that  marine  gravity  results  are  only  as 
good  as  the  cross -coupling  corrections,  horizontal  acceleration 
corrections,  and  Eotvos  corrections.  These  matters  will  not  be 
discussed  in  this  paper.  The  problem  which  we  want  to  examine 
here  is  that  of  effectively  removing  the  "no ise -level"  or  pertur- 
bations in  the  free-air  anomaly,  caused  by  the  bottom  topcgraphy, 
eo  obtain  a residual  gravity  field  which  incorporates  only  the 
effect  of  the  deeper  mass  distributions.  The  situation  is  the 
Inverse  of  the  one  encountered  by  the  paleomagne tic  investigators 
where  they  remove  the  regional  so  that  they  can  look  at  the  local 
effects.  The  most  difficult  topography  effect  to  remove  is  one 
Chat  is  truly  three-dimensional;  one  of  the  most  common  features 
in  marine  bottom  topography  is  the  submerged  seamount.  We 
thought  that  it  would  be  profitable  to  study  the  difference 
between  the  three-dimensional  (3D),  two-dimensional  (2D)  and  one- 
dimensional  (ID)  gravity  effect  of  conical  seamounts.  • 


THE  GRAVITY  EFFECT  OF  A CONICAL  SEAMOUNT 

We  felt  that  it  was  reasonable  to  assume  that  the  form  of  a 
submerged  seamount  could  be  approximated  as  a perfect  cone, 
although  we  realized  that  no  such  seamount  probably  exists.  Our 
chief  purpose  in  addition  to  producing  a table  of  gravity  effects 
was  to  compare  one  and  two-dimensional  approximations  with  the 
three-dimensionally  derived  gravity  effects.  We  examined  the 
bathymetry  character  of  more  than  100  "reasonable"  seamounts  in 
the  central  Pacific  Ocean.  The  slope  angles  were  calculated  from 
the  height  divided  by  base  radius,  and  were  found  to  range  from 
4.5  to  17  degrees  with  10  degrees  as  a reasonable  mean.  Appendix 
A gives  a summary  of  the  developments  leading  to  a group  of 
aquations  used  for  calculating  the  off  axis  vertical  gravity 
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effect  of  a horizontal  circular  disk.  The  seamount  was  then 
approximated  from  stacking  up  to  100  varying  sized  disks.  We 
used  a density  of  2.3  g/cc*,  or  a differential  density  of  (2.3- 
1.027)  g/cc. 

Tables  1,  2 and  3 give  the  vertical  gravltional  attraction 
of  conical  seamounts  of  density  2.3  g/cc  and  having  slopes  of 
5,  10  and  15  degrees  respectively.  The  tables  give  the  radial 
distance  in  km  from  the  axis  corresponding  to  an  even  10  mgal 
unit.  For  cases  where  the  on-axis  mgal  effect  is  not  an  even 
10  mgal  unit,  the  mgal  effect  is  given  in  parentheses  in  place 
of  the  zero-km  number.  The  purpose  in  presenting  the  tables  in 
this  form  is  to  permit  other  investigators  to  construct  graphs 
of  their  own  to  any  scale.  Drawn  as  plan  view  circles  at  the 
appropriate  radii,  they  would  permit  the  gravity  effect  profile 
of  an  off-axis  ship  track  to  be  determined.  If  some  other 
density  p besides  2.3  is  desired, multiply  the  mgal  effect  by 
(p  - 1.027)/(2.3-1.027). 

The  3 -dlmens lonal  results  were  compared  not  only  with 
Talwanl . e t . al . (1959)  2D  equivalents  of  ship  tracks  across  off- 
axis  slices  of  the  cone,  but  also  with  ID  corresponding  approxi- 
mations along  each  track.  The  ship  track  distance  Increments 
were  5 km,  because  this  is  about  the  distance  a ship  will  travel 
in  15  minutes  at  10.5  knots,  and  our  ship  gravity  results  are 
given  for  even  quarter  hour  times.  Table  4 presents  the  com- 
parisons for  the  on-axi3  ship  track  cases.  In  order  to  discuss 
off-axis  cases,  it  seemed  best  to  pick  a few  examples  of  reason- 
able situations  encountered  at  sea.  Figure  1 shows  a half- 
profile of  the  vertical  gravity  effect  of  10'  as  well  as  15' 
conical  seamounts  where  the  base  depth  is  5 km  and  the  depth  to 
the  top  is  one  km.  Figure  2 shows  the  mgal  error  in  using  a ID 


>>  "Figure  1.  Vertical  gravity  effect  at  sea-level  of  two  conical 
seamounts.  Both  have  a base  depth  of  5 km  and  a top  depth  of 
1 km.  Solid  lines  are  for  one  with  a 10'  slope  from  the  hori- 
zontal, and  dashed  lines  are  for  one  with  a 15'  slope.  Corres- 
ponding base  radii  are  indicated  by  the  open  triangles.  Effects 
are  shown  along  ship  tracks  crossing  on-axis,  as  well  as  flank 
crossings  at  5 km  successive  offsets  from  the  axis. 

*In  1971,  R/V  KAMA  KEOKI  with  LaCoste  and  Romberg  gravimeter 
S33  crossed  the  Meeker  Ridge.  Furumoto . et . al . (1971)  had 
established  a reversed  refraction  seismic  line  located  100  km 
east  of  the  axis  of  the  ridge  and  100  km  northeast  of  the  gravity 
profile.  The  km  thickness  and  km/sec  seismic  velocity  values  of 
the  upper  three  layers  of  the  sea  floor  (0.1,  2.3;  0.3,  2.9;  1.8, 
3.6)  along  with  the  Ludwig , Naf e and  Drake  (1970)  corresponding 
densities  of  2.06,  2.21  and  2.33  g/cc  yield  a mean  density  of 
2.30  g/cc.  A Ta 1 wa n 1 (1959)  2D  density  profile  was  used,  and 
2.30  was  confirmed  as  the  best  density  to  use.  This  agrees  with 
conclusions  of  Wool  lard  (1951)  and  to  some  extent  with  those  of 
Strange,  et  al.  (1965). 
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flat-slab  assumption  (upper  graph)  in  correcting  to  the  5 km  base 
level,  as  well  as  the  mgal  error  using  a 2D  Talwanl  method  along 
the  track.  It  is  apparent  that,  on  the  axis,  the  ID  error  is 
enormous,  being  487,  too  large,  while  the  2D  is  much  more  reason- 
able at  16 X too  large.  On  axis,  the  ID  as  well  as  2D  methods 
are  blind  to  the  deficiency  of  mass  laterally.  For  the  10*  sea- 
mount and  ID  assumption,  a conventional  marine  Bouguer  anomaly, 
in  which  the  sea  water  is  replaced  by  rock,  would  over-correct 
the  145  mgal  free -air  anomaly  by  69  mgal  and  appear  as  a large 
negative  spike  of  almost  half  the  amplitude  of  the  original 
anomaly.  Clearly,  then,  the  2D  assumption  is  much  better  because 
in  the  worst  case  (on-axis),  the  145  mgal  free-air  anomaly  would 
be  reduced  to  only-23  mgal,  and  at  10  km  off-axis,  to-8  mgal. 

The  2D-3D  graphs  for  cases  other  than  that  shown  in  Figure  2 are 
quite  similar  in  shape,  and  all  exhibit  the  equality  of  2D  with 
3D  methods  when  the  off-axis  trackline  distance  is  about  60  to 
70  percent  of  the  base  radius.  This  percentage  is  given  in  the 
X RAD  2D  « 3D  column  in  Table  4.  Of  course,  the  2D  assumption 
under-corrects  for  topography  when  the  off-axis  distance  is  too 
. great. 


Figure  2.  Mgal  errors  of  ID  (upper)  and  2D  (lower)  assumptions 
relative  to  the  curves  shown  in  Fig.  1.  Only  the  10*  case  is 
shown  in  the  2D-3D  graph.  Insert  numbers  give  km  offsets  from 
the  axis. 

. EFFECT  OF  FEATURES  OFF  THE  TRACKLIHE 

It  is  important  to  know  th  magnitude  of  the  gravity  effect 
of  features  off  to  the  side  of  the  trackline  and  not  seen 
directly  by  the  echo -s ounding  device  on  the  ship.  The  half- 
angle of  a conventional  echo-sounder  cone  of  energy  is  30* 
relative  to  the  vertical.  This  means  that  the  echo-sounder's 
lateral  effectiveness  ranges  from  3.5  km  in  6 km  of  water  to 
1.7  km  in  3 km  of  water.  From  the  base  radii  given  in  Table  4 
it  is  seen  that  for  all  cases  given  in  Tables  1,  2 and  3,  the 
worst  error  encountered  in  ignoring  seamounts  off  to  the  side  of 
the  track  is  13  mgal  for  a 15“  slope  at  zero  km  to  the  top  and 
6 km  base  depth.  Almost  all  other  cases  are  10  mgal  or  less. 

The  validity  of  ignoring  lateral  topography  would  not  be  true, 
however,  if  there  were  a large  2D  bathymetric  feature  parallel  to 
the  trackline  but  just  out  of  range  the  echo-sounder.  Thus 
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vertical  gravitational  attraction,  at  sea  level  and  on  the  axis,  of  a seamount, 

the  comparison  with  one  as  well  as  two-dimensional  approximations.  Density  is  2.3  g/cc. 
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scraps,  trenches,  ridges,  island'  chains  and  continental  slopes, 
for  example  could  not  be  ignored.  However,  such  situations  would 
be  treated  as  special  cases. 


MOVING  WINDOW  TOPOGRAPHY  CORRECTION 

We  have  developed  a computer  program  which  we  call  TCFAA 
(Topography  Corrected  Free  Air  Anomaly).  Its  purpose  is  to  be 
a relatively  inexpensive  method  of  correcting  the  marine  observed 
free-air  anomalies  for  the  effect  of  bottom  topography  along  the 
track.  A "boxcar  window"  is  centered  on  the  observation  point. 
The  mean  water  depth,  maximum  water  depth,  and  mean  free-air 
anomaly  are  calculated  for  only  the  sequence  of  track  data  within 
the  window.  The  Taiwan!  2D  effect  of  the  topography  relative  to 
the  mean  depth  is  calculated  and  subtracted  from  the  free-air 
anomaly.  Except  for  discrepancies  between  the  2D  assumption  and 
3D  configurations,  the  result  approximates  a "complete"  free  air 
anomaly  relative  to  the  mean  depth  for  the  window.  The  window 
is  advanced  along  the  track  along  with  each  new  observation 
point. 

The  results- of  a test  of  the  program  are  shown  in  Figure  3. 
The  track  was  taken  from  Leg  8 of  KANA  KEOKI  in  1973  and  goes 
from  5 . 5°  S 168°W  (left)  to  3“S  175. 5“W  (right).  All  bathymetric 
"bumps"  are  crossings  of  essentially  conical  seamounts.  The 
center  of  the  track  is  about  120  km  south  of  Canton  Island  and 
the  bathymetric  depression  past  800  km  is  a western  extension 
of  the  Nova  Canton  Trough  and  is  therefore  two-dimensional. 

Figure  3a  shows  original  data.  Figure  3b  shows  the  results 
of  using  a window  width  of  110  km  corresponding  to  1°,  and 
Figure  3c  shows  the  results  of  using  a window  width  of  330  km 
corresponding  to  3"  geographically. 


It 


DISCUSSION  OF  RESULTS 

Figure  3b  shows  a very  close  correspondence  between  the 
mean  free-air  anomaly  and  TCFAA.  This  would  be  expected  because 
TCFAA  is  done  relative  to  the  mean  water  depth.  Over  1°  or  110 
km  the  area  should  be  only  50  percent  isos tatically  compensated, 

. as  has  been  pointed  out  by  Wool lard  (1962) . Because  the  lateral 
v effect  of  seamounts  not  seen  by  the.  echo -sounder  can  be  effec- 

tively ignored,  we  feel  that  the  mean  topography  and  free -air 
anomalies  obtained  with  a HO  km  window  and  a single  ship  track 
should  provide  powerful  information  for  determining  a local  area 
basic  predictor.  This,  in  turn,  can  be  used  to  determine  the 
local  l*xl*  mean  anomaly.  The  TCFAA  profile  for  a 110  km  window 
probably  cannot  be  used  profitably  for  geological  structure 
interpretation.  This  is  because  the  mean  water  depth  does  not 
truly  represent  the  most  commonly  occurring  water  depth  for  the 
whole  region  in  general.  The  most  common  depth  would  be  given 
by  the  "mode,"  which  might  not  be  possible  to  find  in  cases  of  a 
monocline  or  sawtooth  topography,  for  example.  There  is  a way 
out  of  this  problem,  however.  Woollard  (1962)  maintained  that 
topographic  features  must  be  about  3°x3*  or  larger  in  order  to 
be  compensated  more  than  8071.  This  is  the  reason  we  also  show 
results  obtained  with  a 330  km  window. 

Figure  3c  shows  that  the  TCFAA  from  a 330  km  window  should 
be  very  useful  for  local  crustal  structure  interpretations.  In 
particular,  the  departures  of  TCFAA  from  the  330  km  mean  free- 
arca  anomaly  indicate  areas  of  probably  structural  or  tectonic 
interest.  The  mean  free-air  anomaly  as  well  as  mean  topography 
change  very  slowly  and  would  be  quite  useful  for  a regional  basic 
pred ictor . 
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CONCLUSIONS 


1.  Even  in  the  worst  case  where  a ship  track  crosses  the 
axis  of  a submersed  seamount,  the  2D  approximation  for 
removal  of  bottom  topography  gravity  effects  is  consid- 
erably better  than  use  of  a ID  flat-slab  approximation. 

2.  Seamounts  off  to  the  side  of  the  ship  trackline  and  not 
seen  by  the  echo-sounder  can  be  effectively  ignored  as 
regards  the  correction  of  the  free-air  anomaly  for  their 
topographic  effect. 


The  moving  window  TCFAA  computer 
window  can  be  used  very  effective 
to  give  good  l*xl°  free  air  anoma 
addition,  the  program  is  very  use 
crustal  structure  problems  when  u 
dow  because  the  mean  topography  a 
mean  free-air  anomalies  represent 
state"  conditions  while  the  TCFAA 
local  deeper  structural  gravity  e 
contaminated  or  screened  by  the  1 
level." 


program  with  a 110  km 
ly  for  long  ship  tracks 
ly  approximations.  In 
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sed  with  a 330  km  win- 
nd 

regional  or  "steady 

results  represent 
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The  330  km  means  should  be  very  useful  for  setting  up 
regional  basic  predictors. 

The  TCFAA  program  uses  only  what  it  can  see,  and  suffers 
only  from  use  of  a 2D  assumption.  Therefore,  it  is  in 
reasonable  harmony  with  the  principles  of  minimum 
as  tonishment. 
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APPENDIX  A 

THE  VERTICAL  GRAVITATIONAL  ATTRACTION  OF  A 
HORIZONTAL  CIRCULAR  DISK 

The  vertical  gravitational  attraction  of  a very  thin 
horizontal  circular  disk  can  be  computed  from  the  solid  angle 
subtended  by  the  disk.  There  are  several  different  methods  of 
calculating  the  solid  angle,  and  the  computational  efficiency 
of  each  method  vatya  with  the  horizontal  as  well  as  vertical 
distances  from  the  disk. 


The  vertical  gravitational  attraction,  g2,  of 
disk  is  given  by  the  volume  Integral 


8 


z 


Jli£o£!  dv 


horizontal 


(1) 


where  y is  the  gravitational  constant,  6 is  the  density,  r is 
the  distance  of  the  field  point  from  the  source  point  on  the 
disk,  and  6 is  the  angle  between  r and  the  normal  to  the  disk. 
This  is  shown  in  Figure  Al. 


Figure  Al.  Diagram  of  vertical  gravitlonal  effect  of  a hori- 
zontal disk. 


For  h « r,  (1)  reduces  to 

8z  “ v6h  dA  (2) 

A r 


where  A represents  the  area  of  the  disk.  The  solid  angle  at 
field  point  P subtended  by  the  disk  is 


n 


fn  • dA 

J.  r2 


(3) 


wherein  is  the  unit  vector  normal  to  the  surface  of  the  disk 
and  dA  is  the  element  of  area  (dA)  vector  directed  along  the 
line  of  length  r.  Now  (3)  reduces  to 
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which  is  the  same  Integral  as  the  one  In  (2).  Then  (1)  for  a 
thin  disk  can  be  expressed  as 


g^  "•  ybhO  . 

The  value  of  0 on  the  axis  (p“o)  reduces  to 


0 


2n 


(l-2//p-+—7) 


(5) 

(6) 


The  evaluation  of 
more  complicated, 
slons,  and  ochers 
terms  of  Elliptic 


the  Integral  for  a point  off  the  axis  Is  much 
Some  methods  have  used  power  series  expan- 
hive determined  closed  form  solutions  in 
integrals  of  the  first,  second  and  third  kinds. 


Jaf f ey  (1954)  lists  results  of  Zumualt  obtained  by  expan- 
ding the  Integrand  in  a binomial  series  and  then  integrating 
term  by  term.  Since  the  integrand  of  (4),  cos8/r2,  is  a solu- 
tion af  Laplace's  equation,  the  solid  angle  can  also  be 
expressed  in  terms  of  a Legendre  polynomial  series  of  the  form 


n - r An  (~)n+1P  (cos 9) , r>a  (7a) 

nmo 

0 - 2 An  (^Vn  (cose),  r<a  (7b) 

n-n  * 

Where  cos0  - Z/J z2+ ?2 

The  coefficients  An  can  be  solved  for  by  comparison  with  the 
terms  of  the  Zumwalt  expansion  for  the  on-axis  case. 


Refer  to  Jeans  (1958),  p.  431,  for  an  example  of  this 
procedure.  The  particular  result  most  suitable  for  our  purposes 

is 


0 - -2n  Z ?2n  (0)  P2n_i  (cose)  (^-) 

n-1  r 


a.  2n 


r>a 


(8) 

( J-7b) 


which  converges  rapidly  for  r » a.  Other  types  of  series 
expansion  results  are  also  given  by  Jaffey.  He  denote  (8)  as 
our  Equation  (J-7b)  where  prefix  J refers  to  Jaffey. 

Other  approaches  to  the  problem  have  provided  exact  solu- 
tions in  terms  of  the  complete  and  incomplete  Elliptic  integrals. 
One  such  solution  is  derived  by  Paxton  (1959).  Starting  from 
(4)  he  obtains 

0-2 


Refer  to  Figure  A2.  The  line  integral  over  q>  is  obtained  first. 


Figure  A2.  Quantities  involved  in  Paxton's  solution 
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After  expressing  <ps  and  0 In  terms  of  5 and  other  parameters  the 
resulting  Integrals  become  the  complete  elliptic  Integrals  of 
the  first  and  third  kinds.  To  simplify  the  expression,  the 
Elliptic  integral  of  the  third  kind  Is  expressed  In  terms  of 
Heuman's  lambda  function,  A . The  following  results  were 
obtained  by  Paxton:  ° 


0 - 2-rr  - , 23-  K(k)  - vK  (§,k),  p<a  (10a) 

St2  + (p+a)2 

0 - fl  - --  - jZ-  — K(k)  . p-a  (10b) 

/Z2  + (p+a)2 

0 - -====?!===  K(k)  + nA  (5,k)  , p>a  (10c) 

+ (p+a)2  ° 


where 


k2  - 


4pa 


Z2  + (p+a)2 


s in 


-1  ri-k  /«  1 2 

L — 2~ J : 

L 1-k  J 


1/2 


4pa 

r 

(p+a) 


and  K(k)  is  the  complete  Elliptic  Integral  of  the  first  kind. 

The  Heuman  lambda  function  can  be  expressed  in  a series  con- 
taining K(k) , E(k)  (the  complete  Elliptic  integral  of  the 
second  kind),  E , and  k.  To  evaluate  Paxton's  results,  refer  to 
Nagy  (1965)  as  well  as  Abramowltz  and  S tegun  (1970)  , chapter  17. 

Another  exact  solution  for  0 in  terms  of  Elliptic  integrals 
was  obtained  by  Tallqvist  and  referred  to  by  Garrett  (1954). 
Tallqvist  approximated  the  exact  solution  as 

0 - 2n  - 2v  sin9  [l-ik2-i^k4(3-4tan20) . . .]  (11) 

in  this  case  k2  =*  (p2+Z2+a2 -S)  / ( p 2+Z2+a2+S)  , 

2222  221/2  2 222 
S - f(p  +Z  -a  ) + 4a  Z Jj'  sin  0 - (p  +Z2-a  +S)/2S  and 

tan8  ” 2aZ/(S+a2-p  -Z2). 

The  solutions  given  here,  as  well  as  some  others,  were 
examined  to  determine  the  best  series  to  use  in  different 
regions  to  achieve  a given  accuracy  with  a minimum  number  of 
computations  Involved.  The  accuracy  of  (1  was  arbitrarily  set 
atO.OOOl  radian  for  the  test  case.  For  each  solution,  fi  was 
computed  term  by  term  until  the  term  was  less  than  the  desired 
accuracy.  The  resulting  sum  for  each  series  was  then  compared 
with  tabularlzed  solid  angle  values  (flasket,  1957,  1962)  to 
confirm  that  the  series  had  in  fact  converged.  For  computer 
applications,  the  regions  adopted  are  shown  in  Figure  A3  and 
indicate  the  particular  equation  used  and  the  associated  numbers 
of  terms  required  to  achieve  an  accuracy  of  0.0001  radian.  The 
scales  for  p and  Z are  not  linear  as  the  only  purpose  is  to  show 
the  general  regions. 


Figure  A3.  Regions  of  p and  Z adopted  for  0.0001  radian 
accuracy.  Equation  used  and  numbers  of  terms  required  are  given, 
except  that  for  the  Paxton  solution,  14  refers  to  the  number  of 
terms  used  in  Heuman's  lambda  function. 
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